Nature’s Autonomous Oscillators 


Popular Summary 


Nonlinearity is required to produce autonomous oscillations without external time dependent source, 
and an example is the pendulum clock. The escapement mechanism of the clock imparts an impulse 
for each swing direction, which keeps the pendulum oscillating at the resonance frequency. Among 
nature’s observed autonomous oscillators, examples are the quasi-biennial oscillation and bimonthly 
oscillation of the Earth atmosphere, and the 22-year solar oscillation. The oscillations have been 
simulated in numerical models without external time dependent source, and in Section 2 we summarize 
the results. Specifically, we shall discuss the nonlinearities that are involved in generating the 
oscillations, and the processes that produce the periodicities. In biology, insects have flight muscles, 
which function autonomously with wing frequencies that far exceed the animals' neural capacity. 
Stretch-activation of muscle contraction is the mechanism that produces the high frequency oscillation 
of insect flight, discussed in Section 3. The same mechanism is also invoked to explain the 
functioning of the cardiac muscle. In Section 4, we present a tutorial review of the cardio-vascular 
system, heart anatomy, and muscle cell physiology, leading up to Starling’s Law of the Heart, which 
supports our notion that the human heart is also a nonlinear oscillator. In Section 5, we offer a broad 
perspective of the tenuous links between the fluid dynamical oscillators and the human heart 
physiology. 
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ABSTRACT 

Nonlinearity is required to produce autonomous 
oscillations without external time dependent 
source, and an example is the pendulum clock. 
The escapement mechanism of the clock im- 
parts an impulse for each swing direction, which 
keeps the pendulum oscillating at the resonance 
frequency. Among nature’s observed autono- 
mous oscillators, examples are the quasi-bien- 
nial oscillation of the atmosphere and the 22- 
year solar oscillation [1], Numerical models si- 
mulate the oscillations, and we discuss the non- 
linearities that are involved. In biology, insects 
have flight muscles, which function autono- 
mously with wing frequencies that far exceed 
the animals' neural capacity. The human heart 
also functions autonomously, and physiological 
arguments support the picture that the heart is a 
nonlinear oscillator. 

Keywords: Mechanical Clock; Quasi-Biennial 
Oscillation; Bimonthly Oscillation; Solar Dynamo; 
Insect Flight; Human Heart 

1. INTRODUCTION 

Among nature’s oscillators, the lunar tide is a well- 
understood familiar example. The sea level rises and falls 
with a period of about 12 hours, which is produced by 
the gravitational interaction of the moon circling the 
Earth. Another example is the diurnal atmospheric tide, 
which is produced by variable solar heating. Generated 
by external time dependent forcing, the tides are not 
autonomous oscillators. Oscillations that are forced by a 
time dependent source can be understood in terms of a 
linear system. Take n linear equations with n unknowns, 
A.(n, n)-X(n) = S(n), A the square matrix of terms that 
account for dissipation, X the unknown linear column 
vector, and S the source vector. With X(n) = A _1 (n, 
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n)*S(n), A" 1 the inverse matrix, it follows that an oscilla- 
tion X ^ 0 is produced with a source S ^ 0 (Finite dif- 
ferencing produces the matrix elements for A to obtain a 
solution for a system of linear differential equations). In 
space physics, there is a large scientific body of literature 
featuring linear mathematical models that can generate 
and reproduce a variety of observed oscillations applying 
known excitation sources. 

The above linear analysis demonstrates, with X(n) = 
A _1 (n, n)-S(n), that for S = 0, X = 0. This means that 
without an external time dependent source, a linear sys- 
tem cannot produce an oscillation. Autonomous oscilla- 
tors must be nonlinear — and a familiar example of such 
an oscillator is the pendulum clock. Once initiated, the 
clock continues oscillating at the resonance frequency 
that is determined by the length of the pendulum. The 
external source for the clock is the steady force of the 
weight pulling on the drive chain that activates the es- 
capement wheel. With sufficient energy from the sus- 
pended weights to overcome friction, the escapement me- 
chanism provides the impulse/nonlinearity that keeps the 
pendulum oscillating, as illustrated in Figure 1 . 

In mathematical terms, the variables (x, y) interact in a 
nonlinear system, like arxy + bi*xy = Si, a 2 -xy + b 2 -xy = 
s 2 , or for example aj*x 2 + bj*y 3 = si, a 2 *x 2 + b 2 -y 3 = s 2 , 
where a J2 and b i2 account for dissipation such as friction, 
and Si, s 2 are external sources that would be zero for the 
mechanical clock. A simple example of nonlinear inter- 
action is illustrated with the temperature feedback of 
snow. In winter at high latitudes, as it gets colder, the 
precipitation turns into snow. White snowflakes, cover- 
ing the ground, reflect incoming solar radiation to make 
it still colder. The interaction between diminishing radia- 
tion and resulting formation of ice crystals, two variables 
of climate change, say x, y, produce the nonlinear term, 
xy, that accelerates cooling. Another example of climate 
change is the melting of glaciers due to C0 2 -induced 
global warming. As the glaciers melt, the dark surface 
gets exposed and absorbs more solar radiation to ace- 
lerate the warming. Emptying a bottle with fluid would 
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Figure 1 . In the mechanical clock, the escapement mechanism 
produces the impulse/nonlinearity that generates the oscilla- 
tion without external time dependent source. Illustrated are the 
impulse sequence and broadband frequency sequence, which 
produce the pendulum oscillation with resonance frequency. 


be a hands-on experiment demonstrating a nonlinear os- 
cillation. The fluid pouring out, with sufficient speed, is 
compressed in the bottleneck. The resulting increase in 
pressure caused by the flow, a nonlinear feedback, slows 
down the flow to produce a glugging oscillation. 

[1] reviewed the properties of observed fluid dynami- 
cal autonomous oscillators in space physics, foremost the 
quasi-biennial oscillation in the Earth’s atmosphere with 
a period of about two years, and the 22-year oscillation 
of the solar dynamo magnetic field. The oscillations have 
been simulated in numerical models without external 
time dependent source, and in Section 2 we summarize 
the results. Specifically, we shall discuss the nonlineari- 
ties involved in generating the oscillations, and the pro- 
cesses that produce the periodicities. 

Stretch- activation of muscle contraction is the mecha- 
nism that produces the high frequency oscillation of 
autonomous insect flight, briefly discussed in Section 3. 
The same mechanism is also invoked to explain the 
functioning of the cardiac muscle. In Section 4, we pre- 
sent a tutorial review of the cardio-vascular system, heart 
anatomy, and muscle cell physiology, leading up to Star- 
ling’s Law of the Heart, which supports our notion that 
the human heart is also a nonlinear oscillator. In Section 
5, we offer a broad perspective of the tenuous links be- 
tween the modeled fluid dynamical oscillators and the 
human heart physiology. 

2. FLUID-DYNAMICAL OSCILLATORS 

In space physics, like in other areas of environmental 
science, mathematical models have been widely used to 
provide a physical understanding of the observations. 
Mentioned in the introduction, a large number of phe- 
nomena can be understood in terms of linear systems, 
like the diurnal atmospheric tide that is produced by 
variable solar heating. But there is no such external time 
dependent source known, which could produce the ob- 
served quasi-biennial oscillation in the atmosphere or the 
22 -year oscillation of the solar magnetic field. Nume- 
rical models generate the oscillations through nonlinear 


interactions that produce periodicities determined by the 
internal dynamical properties of the fluids. 

2.1. Quasi-Biennial Oscillation (QBO) 

In the zonal circulation of the terrestrial stratosphere at 
low latitudes, the quasi-biennial oscillation (QBO) domi- 
nates, and [2] presented a comprehensive review of the 
observations and physical mechanisms that generate the 
oscillation. Observed with satellites and ground-based 
measurements, the wind velocities vary between east- 
ward and westward regimes that propagate down and 
have periods averaging approximately 28 months. Since 
its discovery more than 40 years ago, research has shown 
that the QBO plays an important role in a multitude of 
processes that affect the climate, as in the transmission of 
solar activity signatures and modulation of atmospheric 
tides for example. In the spirit of the subject matter of 
the present paper, we consider the QBO isolated from 
interactions with other atmospheric oscillations and fo- 
cus on the fundamental processes that generate the QBO. 

In a seminal paper published in 1968, [3] demonstrated 
that the QBO can be generated with planetary waves, and 
their theory was confirmed in numerous modeling stud- 
ies. Subsequent studies showed that the QBO could also 
be reproduced with parameterized gravity waves (GW) 
[4,5], and this is now well accepted [2]. [3] considered 
the QBO being driven by the 6-month semi-annual os- 
cillation (SAO) generated by solar heating. But [6] sub- 
sequently concluded that the SAO was not essential for 
generating the QBO, and this was confirmed in a model- 
ing study [7]. 

From that paper, we present in Figure 2 computer so- 
lutions produced with the Numerical Spectral Model 
[4,8], which incorporates the Doppler Spread Parame- 
terization [9,10] that provides the GW momentum source 
and associated eddy diffusivity/viscosity. Applying factor 
of two different diffusion or dissipation rates, QBO-like 
oscillations are produced with periods of about 33 and 17 
months (Figures 2(a) and (b), respectively), which dif- 
fer in amplitude by about a factor of two. The numerical 
results were generated for perpetual equinox, without 
time dependent solar heating. Without external time de- 
pendent source, the excitation mechanism cannot be lin- 
ear, as demonstrated in Section 1. A nonlinear process 
must be involved in generating the oscillations. 

The nonlinearity is displayed in Figure 3, which shows 
the zonal winds, U, and the accompanying GW momen- 
tum source, MS, varying with altitude. Discussed by 
[1,7], the positive nonlinear feedback associated with 
critical level of wave absorption produces a MS with 
sharp peaks near vertical wind shears, dU/dr. Consider- 
ing the wind oscillation of the QBO, varying with fre- 
quency co , the nonlinear part of the MS, to first order, can 
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Figure 2. Altitude contour plots of equatorial zonal winds, 
positive eastward and negative westward. The winds were de- 
rived with diffusion/dissipation rates a factor of 2 larger in (b) 
than in (a). The resulting oscillation periods and amplitudes de- 
crease with increasing dissipation rates; they are a factor of 2 
smaller in (b) than in (a). A constant gravity wave (GW) flux 
provided the energy for the momentum source that generated 
the QBO-like oscillations. Figure taken from [7]. 

be written in the approximate form [dU(co,r)/dr] 3 , of odd 
power. With complex notation, a source in the form 
[Exp(icjt)] 3 produces the term with Exp[i(a> + co + cu)t] 
for the higher order frequency, 3 co. Such a nonlinear 
source also generates Exp[i(ou + co - co) t] for the funda- 
mental frequency co, which maintains the oscillation. 

2.2. Bimonthly Oscillation (BMO) 

Zonal-mean oscillations with periods around 2 months 
have been seen in ground based and satellite wind 
measurements in the Earth atmosphere [12-14]. The 
Numerical Spectral Model generates such oscillations 
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Arbitrary Units 

Figure 3. Snapshots are shown of zonal winds, and effective 
GW momentum source that generates the oscillations in Figure 
2(b). The peak momentum source occurs near the maximum 
vertical velocity gradient; it varies in phase with the gradient 
and is strongly nonlinear. Figure taken from [11]. 

[15], and numerical experiments show that they are not 
produced when the gravity waves (GW) propagating 
north/south are turned off. In Figure 4(a), a computer 
solution is presented from that paper, which is produced 
only with the meridional (north/south) GW momentum 
source. Regular bimonthly oscillations (BMO) are gene- 
rated with a period of about 2 months, which resemble 
the quasi-biennial oscillation (QBO) above discussed. 

The similarity between the BMO and QBO is also 
evident in Figure 4(b), where meridional winds, V, are 
presented together with the associated momentum source. 
Like the GW source of the QBO in Figure 3, the mo- 
mentum source in Figure 4(b) is sharply peaked near 
meridional wind shears. Discussed in Section 2,1, this 
impulsive momentum source can be viewed approxi- 
mately as a nonlinear function of the velocity field, which 
has the character, [dV(ft>,r)/dr] 3 . Such a non-linearity can 
generate an oscillation without external time-dependent 
forcing. Like the QBO, the meridional BMO can be un- 
derstood as a nonlinear oscillator. 

Like the QBO, the BMO is dissipated by viscosity. But 
unlike the zonal winds of the QBO, the meridional winds 
of the BMO produce pressure variations that counteract, 
and dampen, the winds to produce a shorter dissipative 
time constant. This additional thermodynamic feedback 
explains why the period of the BMO is much shorter 
than that of the QBO. In both cases, the oscillation 
periods are determined by the dissipation rates. 

2.3. Solar Dynamo 

The solar magnetic field is observed varying with a 
period of about 22 years. During periods of enhanced 
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Figure 4. (a) Meridional wind oscillations with periods of 
about 2 months, (b) Snapshots of normalized winds with mo- 
mentum source that produces sharp peaks near vertical wind 
shears, similar to Figure 3. Figure taken from [15]. 


magnetic fields, irrespective of polarity, sunspots form 
around the equator that are associated with enhanced 
solar radiation. The resulting 11 -year cycle of solar ac- 
tivity produces variations that strongly depend on the 
wavelength of the emitted radiation. At lower altitudes in 
the atmosphere near the ground, the solar cycle effect is 
relatively small. But with increasing altitude, the effect 
increases due to the shorter wavelength ultraviolet radia- 
tion that is absorbed, and above 150 km, the medium of 
artificial satellites, the radiative input can vary by as much 
as a factor of 2 during the solar cycle. The solar cycle 
and its periodicity are variable, and empirical models that 
employ the observed precursor polar magnetic field dur- 
ing the minimum have been very successful in predicting 
the magnitude and duration of the solar maximum [16]. 

Pars pro toto, we refer here to the magneto hydro-dyna- 
mic (MHD) model of [17]. [17] generate the solar os- 
cillation with a zonal-mean magnetic field model, which 
produces the solar oscillation shown in Figure 5. In this 
model, the poloidal (radial) magnetic field is wound up 
by solar differential rotation to produce a toroidal (long- 
itudinal) field in the tachocline below the convective en- 
velope. The buoyant toroidal field emerges at the top of 
the convection region to form sunspots. In that process, 
the Coriolis force is twisting the toroidal field to rege- 
nerate the poloidal field. 

[17] artificially add in their equation for the poloidal 
magnetic field a nonlinear term presented here in a Tay- 
lor series expansion with variable where B $ re- 

presents the toroidal magnetic field, (r,&) are spherical 
polar coordinates (r c , for the tachocline region), and B 0 is 
taken to be constant, 

S(r,9;B ,) °c [sin0cos0] [b, -B}/b 2 0 +B 5 Jb 4 0 +■■■] 

( 1 ) 



Figure 5. Shown is the radial magnetic field varying with a 
period of 20 years. Solid (dashed) contours correspond to posi- 
tive (negative) magnetic fields, which are spaced logarithmi- 
cally, with three contours covering a decade of field strength. 
Figure is taken from the MHD dynamo model of [17]. 

The nonlinear source terms of Eq.l have the property 
that they are of odd (e.g., 3 rd ) power, a similar nonline- 
arity also appears in the classical dynamo model of [18]. 
In dynamo models, the nonlinear terms have the char- 
acter of the momentum source (Figure 3) that generates 
the quasi-biennial oscillation (QBO). The QBO and solar 
dynamo have furthermore in common that the oscillation 
periods are determined by the dissipation rates. For the 
QBO, it is the eddy viscosity; and in the MHD model of 
[17], the magnitude of the dissipating meridional circu- 
lation mainly controls the period of the solar cycle. 

3. INSECT FLIGHT 

Insect flight is an outstanding example of autonomous 
oscillation, and the muscle physiology involved has 
served as a model for understanding the human heart. 

In an advanced class of insects, wing frequencies as 
large as 1000 Hz are produced, which are determined by 
the inertia of the wings and far exceed the frequency 
capacity of the animals’ neural system [19]. This auto- 
nomous flight mechanism is produced by stretch- 
activation of muscle contraction [20,21], which is also 
important for the mammalian heart. Insect flight is very 
demanding energetically, and stretch-activation proves to 
be well suited to match the wing frequency to the ae- 
rodynamic load, making the insect “one of the most 
successful animals that have conquered air” [22]. In this 
fascinating class of insects, the flight muscles are inte- 
grated into the animal’s elastic thorax (Figure 6), chang- 
ing its shape. They are therefore referred to as indirect 
flight muscles. The pair of dorsoventral muscles (Figure 
6(a)) contract and pull the tergal plate at the back of the 
insect down, which causes the hinged wings to flip 
upwards. The longitudinal muscles (Figure 6(b)) run 
horizontally and compress the thorax from front to back, 
which causes the thorax to bow upwards and make the 
wings flip downwards. These antagonistic pairs of flight 
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Autonomous Insect Flight 




Figure 6. Insect cross-sections illustrate (a) con- 
traction by the dorsoventral muscles that pull 
the tergal plate down to make the hinged wings 
move up, and (b) contraction of the longitu- 
dinal muscles that push the tergal plate up and 
the wings down. Figures taken from [23]. 

muscles work in tandem through stretch-activation. When 
the dorsoventral muscle contracts to lift the wings up 
(Figure 6(a)), the longitudinal muscle is stretched to 
activate its contraction that pulls the wings down (Figure 
6(b)), and vice versa. The two antagonistic muscles alter- 
nately stretch and contract. Stretch activates the muscle 
contractions that force the wings up and down to produce 
the high frequency oscillation of autonomous insect 
fright. Nerve impulses command the insect flight but at a 
rate much slower than the wing oscillation; hence it is 
referred to as asynchronous flight. 

4. HUMAN HEART 

A tutorial review is presented of the human heart func- 
tion, which provides the framework for the notion that 
the autonomous heart can be understood as a nonlinear 
oscillator. 

4.1. Heart Physiology 

In the cardiovascular system [24], illustrated in Figure 
7, the human heart exchanges blood with the pulmonary 
and systemic circuits that complement each other. Oxy- 
gen-rich and C0 2 -poor blood, returning from the lungs in 
veins, is pumped in arteries by the left side of the heart to 
the tissue cells of the human body. The returning blood 
in veins, oxygen-poor and C0 2 -rich, is pumped in arter- 
ies by the right heart into the lungs that pick up oxygen 
and shed C0 2 . 

Figure 8 illustrates the human heart anatomy. The 
heart functions through the interplay between the receiv- 
ing chambers, the atria, and the main pumping chambers, 


Cardiovascular Circuit 

Right Side Leftside 



Figure 7. Pulmonary (top) and sys- 
temic (bottom) circuits of the cardio- 
vascular system. The left side of the 
heart receives in veins oxygen-rich and 
C0 2 -poor blood (red) from the lungs, 
which is pumped in arteries to the tis- 
sue cells of the body. The right side of 
the heart receives in veins oxygen- 
poor and C0 2 -rich blood (blue) from 
the body, which is pumped in arteries 
to the lungs. 

the ventricles, which are activated by the myocardium 
heart muscles. Oxygen-rich blood from the lungs returns 
in veins, fills the left atrium, and enters through the mi- 
tral valve into the left ventricle. And the left ventricle 
pumps the blood through the aortic valve and arteries to 
the main body tissues. On the right side of the heart, the 
returning oxygen-poor blood in veins fills the atrium and 
enters through the tricuspid valve into the right ventricle, 
which pumps the blood through the pulmonary valve into 
the arteries of the lungs. During the diastolic phase of the 
heartbeat, blood enters the expanding ventricles. During 
the systolic phase, the contracting ventricles then pump 
the blood with increased pressure into the lungs and the 
rest of the body — and Figure 9 describes the physical 
signatures that characterize the cardiac cycle associated 
with the actions inside the left side of the heart. 

As illustrated in Figure 9, the closure of the aortic 
valve, audible in the phonogram, signals the start of the 
diastolic phase of the heart rhythm. Blood from the pre- 
vious cardiac cycle is prevented from returning to the 
heart, and the muscle of the left ventricle is relaxed, oc- 
cupying a minimum volume of about 50 ml (The left 
atrium is also relaxed at reduced volume). With the open- 
ing of the mitral valve, the main diastolic phase of the 
ventricle is activated. Blood returning from the lungs ra- 
pidly flows through the atrium into the left ventricle 
causing it to expand. At some point during this phase of 
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Figure 8. In the human heart, the left and right atria are filled 
with blood returning in veins (equipped with valves) from the 
lungs and main body tissues, respectively. Blood is pumped 
through the mitral and tricuspid valves into the corresponding 
ventricles. The ventricles pump the blood through the aortic 
and pulmonary valves into the main body tissues and lungs, 
respectively. The pumping actions are provided by the myo- 
cardium muscles surrounding the atria and ventricles. Figure 
taken from Wikipedia. 


expansion, the left atrium begins to contract, and the 


atrium pressure increases. During the atrial systole, blood 
is actively pumped through the mitral valve into the ven- 
tricle, causing it to expand further. Eventually, the pres- 
sure in the ventricle becomes large enough to counter the 
atrium pressure, and this brings the incoming flow to a 
halt, causing the mitral valve to close. Under enhanced 
pressure, the ventricle has expanded to the peak volume 
of about 130 ml, almost three times the value at the be- 
ginning of the diastolic phase. 

As shown in Figure 9, the systolic phase of the ventri- 
cle starts with the closing of the mitral valve, audible in 
the phonogram. At this point, the blood pressure abruptly 
increases as the ventricle contracts, evident in the de- 
creasing volume. Following that event after a short pe- 
riod of time, the pressure is large enough to open the 
aortic valve, and the blood is ejected into the human 
body. At the end of the systolic phase, the ventricle then 
relaxes with decreasing volume, and the blood pressure 
goes down. Below a certain value of the ventricle pres- 
sure, the aortic pressure closes the aortic valve, and the 
systolic phase comes to an end. 

The above-described features of the cardiac cycle ap- 
ply to the left heart, which supplies the human body with 
oxygen-rich blood received from the lungs. The same 
kinds of features are also involved in the functioning of 
the right heart, which delivers to the lungs the oxy- 
gen-poor blood taken from the rest of the human body. 
The physical demands on the pumping actions of the 
right atrium and ventricle though are weaker than those 
on the left side. 


Rapid inflow Ejection 



Figure 9. For the left side of the heart, the time variations for the observed pressure and 
volume of aorta, atrium and ventricle. Also shown are the electrocardiogram and phonogram 
signals. Different phases of the cardiac cycle are marked by the opening and closing signals of 
the aortic and mitral valves and by the abrupt variations of blood pressure/volume. Figure taken 
from Wikipedia. 
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4.2. Cardiac Pacemaker 

In contrast to the skeletal muscle, the cardiac muscle 
has its own conduction system and is not controlled by 
central nerves. The autonomic portion of the peripheral 
nervous system can control the human heart, such that 
sympathetic nerves initiate the cardiac response to “flight 
or fight”, and parasympathic nerves calm it down. Basi- 
cally, the human heart functions autonomously, and its 
auto rythmic contractions occur spontaneously. Located 
in the right atrium, the muscle cells of the sinoatrial (SA) 
node act like a cardiac pacemaker. The electric impulses 
generated by the SA node are associated with the con- 
traction of the atria. After a short delay, the impulses then 
spread through the atrioventricular node (AV), located 
between the right and left atria, and then reach a large 
group of muscle fibers enclosing the right and left ven- 
tricles. These impulses are associated with the contrac- 
tion of the ventricles. The impulse transmission through 
the heart’s conduction system is detected in electrocar- 
diograms (Figure 9), where the P wave signals the atrial 
contraction, and the QRS wave signals the delayed con- 
traction of the ventricle. 

The cardiac pacemaker displays features that resemble 
the nonlinear clock mechanism illustrated in Figure 1. In 
the mechanical clock, the escapement mechanism pro- 
duces an impulse for each swing direction to keep the 
pendulum oscillating at the resonance frequency. Like 
the nonlinear features of the clock mechanism, nonlin- 
earity is evident in the impulse sequence of the EKG and 
in the abrupt variations of the ventricle blood pressure/ 
volume (Figure 9), discussed in the literature [25-32]. 
The question is what is the physiological mechanism that 
produces the nonlinear signatures of the autorythmic 
heart function. 

4.3. Muscle Cell Physiology 

The sarcomere is the basic building block of the mus- 
cle cell, and its physiology has recently been reviewed 
[33]. Illustrated in Figures 10(a) and 10(d), the principal 
protein structures of the sarcomere are the elongated 
thick myosin and thin actin protein filaments. First ob- 
served with high resolution electron microscopes [34,35], 
the two filaments in the stretched and extended sar- 
comere (Figure 10(a)) slide past each other to produce 
the shortened configuration (Figure 10(d)) of the muscle 
contraction. In this sliding action, the titin protein acts 
like a sponge. 

At the molecular level, the mechanism of muscle con- 
traction is illustrated in Figures 10(b) and (c). During 
contraction, the myosin heads attach to the actin proteins 
and pull them towards the middle of the sarcomere. Actin 
filaments slide past myosin filaments, and the sarcomere 
shortens. The process is regulated by two proteins, tro- 
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Figure 10. The sarcomere is made up of elongated actin and 
myosin protein filaments, shown for the stretched/relaxed (a) 
and contracted (d) muscle. In the stretched sarcomere, the tro- 
pomyosin protein covers the actin and blocks the myosin in- 
teraction (b). Ca^ and troponin protein combine to dislodge the 
tropomyosin. This opens the actin globules to the myosin inter- 
action, which is produced by cross-bridges that attach to the 
actin filament and drag it to the center of the sarcomere, hence 
shortening it (c). Adenosin triphosphate (ATP) provides the en- 
ergy for the sarcomere contraction. Figures taken from Wiki- 
pedia, and www.ucl.ac.uk. 


pomyosin and troponin. Tropomyosins are coiled pro- 
teins that run along the actin filament, where they can 
cover the binding sites for myosin and thus block the 
muscle contraction (Figure 10(b)). Calcium from the 
sarcoplasma is required for the muscle contraction. The 
Ca^ ions bind with the troponin complex, which recon- 
figures the tropomyosin such that the globular actin pro- 
teins become exposed to the myosin heads (Figure 
10(c)). Myosin heads then deliver the mechanical force 
for muscle contraction. That force is produced by so 
called cross-bridges [36,37], which are hinged elements 
of the myosin heads that bend and attach to actin, and 
then drag the actin filament to produce muscle contrac- 
tion. The energy required for this force is provided by 
adenosin triphosphate (ATP) molecules that attach to 
myosin heads. ATP is unstable in water. When ATP is 
hydrolized, the bond connecting the phosphate is broken 
to produce adenosin diphosphate (ADP), and energy is 
released. The chemical energy thus liberated powers the 
bending myosin heads that push the actin filaments (Fig- 
ure 10(c)) to produce the sarcomere muscle contraction. 
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4.4. Starling’s Law of the Heart and 
Stretch-Activation of Muscle 
Contraction 

We believe that the nonlinearity of the heart function 
is manifest in the stretch-activation of cardiac muscle 
contraction, consistent with Starling’s Law of the Heart. 

In a series of laboratory experiments with dogs led by 
Starling [38], the team of scientists observed that the 
increasing blood in veins filling the cardiac ventricle led 
to increased stroke volume caused by muscle contraction. 
This property of the cardiovascular system is found in 
finimals and humans, and it is independent of neural in- 
fluences on the heart. In mechanical terms. Starling’s law 
implies that the initial stretching of the cardiac muscle 
affects contraction. The more the muscle is stretched the 
greater the force of contraction. As a result, the heart 
naturally adapts to the physiological condition and de- 
mands of the body, Starling’s Wisdom of the Body. 
Quoting Wikipedia, Starling’s Law of the Heart states 
that the stroke volume of the heart increases in response 
to an increase in the volume of the blood filling the heart 
( the end diastolic volume). The increased volume of the 
blood stretches the ventricular wall, causing the cardiac 
muscle to contract more forcefully. 

Stretch-activation of muscle contraction is a pronoun- 
ced, and outstanding, feature of insect flight muscles [20, 
21], and significant progress has been made explaining 
the mechanism at the molecular level, reviewed by [22] 
and briefly illustrated in Section 4.3. Applying electron 
microscopy, [39] studied the force production of myosin 
cross-bridges (Figure 10(c)), separating the effects of 
Ca^ 4 abundance and stretching of the sarcomere. For this 
experiment, the giant water bug Lethocerus was used, 
where isolated muscle fibers can oscillate for several 
hours. The experiment clearly demonstrated that stretch- 
activation and Ca 4 ^ ions are complementary pathways for 
cross-bridge attachment and force production. The data 
suggested that stretch-induced distortion of attached 
cross-bridges relieves the tropomyosin blocking (Figure 
10(b)) of binding sites on actin, so that the maximum 
force of contraction can be generated at low calcium con- 
centrations. Following up on this study, [40] discovered 
new cross-bridges between myosin and actin filaments, 
involving troponin, they named troponin bridges. They 
exposed the oscillating muscle to intense X-rays and re- 
corded during muscle contraction the diffraction patterns 
associated with the positions of myosin heads, actin fila- 
ments, troponin, and tropomyosin. Based on these data, 
they proposed a comprehensive, self-consistent structural 
model for stretch-activation. When the muscle fibers are 
stretched, provided there is some Ca 44 present, the tro- 
ponin bridges assist in pulling tropomyosin away to ex- 
pose the myosin binding sites on actin, thus activating 
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contraction (Figure 10(c)). While the thin actin filaments 
are opened due to stretch, the thick myosin filaments are 
twisting in response to stretch. The result is that more 
myosin heads are brought to actin target zones during 
muscle stretch so that more cross-bridges are formed for 
larger force production. 

In the skeletal muscle, nerve impulses stimulate the re- 
lease of Ca 44 , which prompts the muscle contraction. 
Calcium is also important for the cardiac muscle but does 
not control its operation. Discussed by [41], the rhythmic 
nature of the cardiac muscle contraction calls for a func- 
tional analogy with stretch-activation of the insect flight 
muscle. Moreover, stretch-activation is called for by the 
steep length tension relationship of the cardiac muscle, 
which is necessary to explain Starling’s Law of the Heart 
[42,43], In cardiac muscles, the importance of stretch- 
activation has been demonstrated [44-46]. [47] investi- 
gated the stretch-activation of cardiac muscle contraction 
and its dependence on the Ca 44 concentration. Calcium 
levels were varied, and the force response to stretch was 
evaluated. The data showed that stretch-activation is 
most pronounced at low levels of calcium activation 
when the thin actin filament offers many sites for the 
formation of additional cross-bidges that can produce 
muscle contraction. They concluded that stretch-activa- 
tion produces strong-binding cross-bridges, which in turn 
open up target zones for new cross-bridges on the actin 
filament, and thus increase the force response during 
muscle contraction. Another force producing pathway 
was pointed out by [48]. They proposed that the Ca 24 
dependent sensitivity of muscle contraction can be mo- 
dulated by the stretch or length dependent variations in 
the lateral spacing between actin and myosin filaments, 
which affects the formation of myosin cross-bridges. At 
the molecular level, [49] employed X-ray diffraction to 
determine the cross-bridge proximity of the myosin heads 
in relation to the thin actin filament. It showed convinc- 
ing evidence for a reduction of spacing between the thick 
and thin filaments, which increased the number of cross - 
bridges involved in the force production during the early 
phase of the ventricular contraction. As shown for the 
autonomous operation of the insect flight muscles, there 
is now abundant observational evidence that stretch-acti- 
vation is of central importance for the cardiac muscle 
contraction, consistent with Starling’s Law of the Heart 
[38], 

4.5. Nonlinear Heart Function 

Given that stretch-activation explains the function of 
the cardiac muscle, consistent with Starling’s Law of the 
Heart, what justifies our picture that the heart is a non- 
linear oscillator? 

Referring to the introduction (Section 1), simple ex- 
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amples of nonlinearity are evident in the feedback of 
snow, which reflects the solar radiation to accelerate 
cooling, and in melting glaciers due to C0 2 -induced 
global warming, which further accelerates warming. Non- 
linearity is also apparent in the oscillation of fluid 
streaming out of a bottleneck, being chocked by pressure 
feedback. From the numerical results discussed in Sec- 
tions 2. 1 and 2.2, we learn that the nonlinear interactions 
between gravity waves, x, and winds, y, produce the 
steep nonlinear momentum source, xy, shown in Figure 
3 and Figure 4(b), respectively. 

In the case of the human heart, the diastolic built up of 
the ventricle blood pressure and resulting stretching of 
the muscle sarcomere (Figure 10(a)) could be viewed as 
variable, x. The stretch induced formation of myosin 
cross-bridges (Figure 10(c)), discussed in Section 4.3, 
could be the other variable y, which produces the muscle 
contraction and built-up of systolic blood pressure. Since 
stretch induces muscle contraction, more stretch pro- 
duces more muscle contraction, we believe that the non- 
linear nature of the heart function seems well established. 
Nonlinear signatures of this interaction are seen in abrupt 
variations of the ventricular blood pressure and volume 
(Figure 9), and in the electrocardiogram that signals dif- 
ferent phases of the cardiac cycle (Figure 9), as dis- 
cussed in the literature [25-32]. In our picture of the 
autonomous human heart, the rhythmic oscillation is 
produced by the nonlinearity built into the stretch acti- 
vated contraction of the cardiac muscles, consistent with 
Starling’s Law of the Heart. This physiological mecha- 
nism could prompt the pacemaker signals, seen in elec- 
trocardiograms, which are transmitted through the car- 
diac conduction system to coordinate and synchronize 
the muscle functions. Like the nonlinear clock that needs 
to be activated initially to start oscillating, the human 
heart needs to be pumped up to start functioning again. 

The above stretch-activation of muscle contraction ap- 
plies to the ventricles, the main pumping chambers of the 
heart. The atria that pump the blood into the ventricles 
during the atrial systol, as seen in the blood pressure dia- 
gram (Figure 9), are the complementary antagonists to 
the ventricles. As the muscles in the atria contract, the 
muscles in the ventricles stretch, analogous to the earlier 
discussed stretching and contraction of the antagonistic 
autonomous insect muscles (Section 3, Figure 6). The 
atria could be seen as being stretched by the blood re- 
turning in the venous system of the pulmonary and sys- 
temic circuits of the human body (Figures 7 and 8). And 
stretching of the atria could induce atrial contraction to 
produce the other nonlinear component of the oscillating 
heart, thus complementing the nonlinear stretch activa- 
tion of the ventricles. 

For the heart to oscillate, an antagonistic and comple- 
mentary counter-part is required to oppose the pumping 
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action during the systolic phase. The blood in veins with 
valves (Figure 8), returning from the lungs and the rest 
of the body, must be pumped into the atria during the 
diastolic phase (Figure 9). This pumping action is pro- 
vided by: 1) the contraction of skeleton muscles and 
pulmonary pump that squeeze the elastic veins; 2) the 
contraction of smooth muscles surrounding the enervated 
veins, which is prompted by the sympathetic nervous 
system that responds to the stress and demands of the 
body; and 3) perhaps, the stretch induced contraction of 
the veins’ smooth muscles in response to the built-up of 
blood pressure. 

5. CONCLUDING REMARKS 

The fluid dynamical oscillators discussed have been 
simulated with nonlinear numerical models. In the quasi- 
biennial oscillation (QBO) and bimonthly oscillation 
(BMO), the gravity wave interactions with the atmos- 
pheric circulation produce momentum sources that are 
strongly nonlinear as shown respectively in Figure 3 and 
Figure 4(b). For the solar dynamo model, a nonlinear 
source term (Eq.l) is employed to generate the 22-year 
oscillation. These nonlinear autonomous oscillators were 
not controlled, and constrained, by external time de- 
pendent forcing. The periods of their oscillations are 
solely determined by the internal dynamical properties of 
the fluids. 

Numerical results for the QBO (Figure 2) show that 
the amplitude and period of the oscillation are reduced 
by about a factor of two when the dissipating viscosity is 
increased by that factor. The period of the BMO (Figure 
4) is generated by the thermodynamic feedback of the 
meridional winds. In the solar dynamo model [17], the 
oscillation period of the magnetic field (Figure 5) varies 
inversely with the magnitude of the dissipating meri- 
dional circulation. The physical properties of these fluid 
dynamical oscillators resemble those of the nonlinear 
human heart, which is affected by the constraints and 
demands of the body. With arterio sclerosis due to nar- 
rowing of the arteries for example, the heart experiences 
more resistance in delivering the blood. As a result, the 
heart rate increases, and the volume of blood delivered in 
each stroke decreases — analogous to the changing QBO 
in response to the increased rate of dissipation. 

The nonlinear fluid dynamical oscillators discussed 
were chosen to function autonomously without external 
time dependent source. The mechanical clock was not 
shaken, and the QBO was not forced externally by time 
dependent solar heating. In reality, numerical experi- 
ments show that, under the name of a pacemaker, the 
seasonal variations amplify the QBO and modulate its 
periodicity [50,51]. The human heart does not function 
solely through stretch-activation. The heart functions in 
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concert with the vascular system, where sympathetic 
nerves activate the smooth muscles around enervated 
veins that pump the blood to the heart, in response to 
stress or exercise for example. This autonomous, and 
harmonic, function of the human heart can be interrupted 
by sympathetic nerves that can initiate the cardiac re- 
sponse to “flight or fight”, and by parasympathetic nerves 
that calm down the heart. 

The basic physical processes involved in describing 
the present fluid dynamical oscillators are well under- 
stood. To obtain global-scale numerical computer solu- 
tions, a parameterization algorithm was employed to 
simulate the gravity wave interactions that generate the 
QBO and BMO, and a nonlinear source term was intro- 
duced to generate the solar dynamo oscillation. 

Compared to the nonlinear oscillators in space physics, 
the multifaceted operation of the human heart and 
cardio -vascular system is far more complicated. To deal 
with that complexity, one could envision parameterize- 
tion schemes that describe the observed properties of 
stretch activated muscle contraction at the molecular 
level. Nonlinearity would come into play by allowing 
stretching to influence contraction in the atrial and ven- 
tricle heart chambers (Figure 8). To produce the oscilla- 
tion would require that the heart is made to interact with 
the vascular system (Figure 7). This antagonistic and 
complementary process to the heart function would re- 
quire a parameterization that describes the volume of 
returning venous blood in relation to the stress and de- 
mands of the human body, satisfying Starling’s Law of 
the Heart and Wisdom of the Body. Dissipation would 
have to be built into the model, allowing for example for 
the resistance the blood experiences in the arteries, 
analogous to the viscous dissipation of the QBO and 
friction in the mechanical clock. 

The complex cardio-vascular system could perhaps be 
compared with the weather systems near the ground and 
in space, where a multitude of nonlinear interactions pro- 
duce a complex phenomenology. Physics based numeri- 
cal mathematical models, aided by parameterizations, 
have greatly advanced our understanding of the envi- 
ronment to produce powerful predictive capabilities. Con- 
sidering the enormous advances made in describing and 
understanding the heart physiology, we believe that the 
next decade will probably see mathematical models that 
can simulate the operation of the human heart. 
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